NUCLEAR HYDRODYNAMICS IN THE INNER CRUST 
OF NEUTRON STARS* 

PiOTR Magierski^ and Aurel Bulgac^ 

^Faculty of Physics, Warsaw University of Technology 
ul. Koszykowa 75, 00-662 Warsaw, Poland 
Piotr.Magierski@olimp.if.pw.edu.pl 
^Department of Physics, University of Washington, 
Seattle, WA 98195-1560, USA 
bulgac@phys . Washington . edu 

(Received October 17, 2003) 

In the inner crust of a neutron star, due to the high density and pres- 
sure, nuclei which are still present, are immersed in a neutron superfluid. 
One then expects that the dynamical properties of nuclei are significantly 
affected. In order to estimate the magnitude of the effect associated with 
the presence of a superfluid medium, we formulate the hydrodynamical ap- 
proach to the nuclear dynamics in the inner crust of neutron stars. We 
calculate the renormalized nuclear mass and the strength of the medium- 
induced interaction between nuclei. We argue that these effects noticeably 
modify the properties of the Coulomb crystal in the inner crust. 

PACS numbers: 21.65.-hf, 97.60.Jd, 26.60.-hc, 95.30.Lz 



1. Introduction 

The inner crust of neutron stars has a complex structure which depends 
on the nucleon density. In the outer parts of the inner crust spherical nuclei 
are immersed in a neutron superfluid. They form a crystal lattice stabilized 
by the Coulomb interaction. At larger densities the phase transition to the 
deformed nuclear phase occurs, leading eventually to the predicted existence 
of a great variety of exotic nuclear structures. The static properties of the 
nuclear matter under such extreme conditions have been a subject of a 
considerable theoretical effort [111111311101117111111110111111111131111113 
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In particular, the phase transition pattern at the bottom of the crust 
has been analyzed using methods based on the density functional theory. 

However, in order to analyze the thermal and electric conductivities of 
the crust, which are important for understanding of the thermal evolution 
of neutron stars one has to go beyond the static approximation. Since 
the part of the inner crust is filled with nuclei likely forming a Coulomb 
crystal, the dynamical properties of such systems have been studied by many 
authors (see e.g. |17| and references therein). An important ingredient of 
the theory of a Coulomb crystal is the plasma frequency, which defines the 
energy scale for the lattice vibrations. It depends on the ion density, ion 
charge and masses. Although the two first quantities can be calculated using 
the usual density functional approach, the last one requires some caution, 
since the nuclei are immersed in a fermionic environment and their bare 
masses are renormalized, in order to take into account the interaction with 
the surrounding neutrons. Namely, as nuclei oscillate around the stable 
positions they bring the fermionic medium into motion as well. 

The example we are considering here is a particular realization of the 
problem of a heavy particle (impurity) dynamics in a Fermi system. In 
general, one may expect that the equation of motion for such an impurity 
moving through the fermionic medium could be expressed in the form: 

where g is a coordinate describing the position of the nucleus, r] is the 
dissipation coefficient and Eint denotes the reversible part of interaction 
energy with surrounding fermions. There are two problems associated with 
the above equation. First, one has to determine the mass M^*^" of the 
impurity immersed in the fermionic medium. Certainly a bare mass of an 
impurity should be renormalized, since the impurity travels in the fermionic 
medium in a cloud of medium excitations. The second problem is associated 
with the fact that the simple Newtonian equation of motion with friction 
cannot be translated directly to quantum mechanics. There is no analogue of 
the concept of friction in quantum mechanics as friction is not time-reversal 
invariant. The solution to this puzzle is the coupling of the considered 
impurity to an environment, which consist of various excitations of the Fermi 
liquid. However deriving from this approach an effective equation of motion 
for the impurity is usually based on some serious approximations. There are 
a number of papers on the problem [UlIiaiinilllll^llSlillliliaiZailHl, 
which however have little overlap with the present approach. The reason 
is threefold, first, we consider the rather particular case of an impurity 
(nucleus) which is not only heavy, but also large (i.e. kpR ^ 1, where kp 
is the Fermi momentum and i? is a nuclear radius). Second, the neutron 
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bath in which the nucleus travels is incompressible, and last but not least, 
the system is superfluid, since the temperatures we are interested in, are far 
below the critical temperature for loss of superfluidity, which is believed to 
be of the order of 1 MeV. 

Hence in the following we will make certain assumptions dictated by 
the special properties of the nuclear system. Namely, we assume that the 
neutrons outside nuclei form an incompressible and irrotational liquid. It 
implies that their dynamics can be described by a velocity potential. Second, 
we neglect the energy dissipation effects, assuming that the energy trans- 
fer between the nucleus and the medium is a fully time-reversible process. 
This assumption seems to be well justified for sufficiently low temperatures. 
Although it is clear that at finite temperatures there always exists a gas 
of quasiparticle excitations, which are responsible for energy dissipation, 
however if wc limit our analysis to small temperatures, as compared to the 
critical temperature of the system, we may disregard this effect, since the 
number of quasiparticle excitations decreases exponentially with decreasing 
temperature. Last, but not least, we will assume that the flow velocity in- 
duced by a moving nuclear impurity is below the critical velocity for the 
loss of superfluidity. Again it is justified for low energy lattice excitations 
which are the most interesting in the context of the neutron star crust. 

2. Hydrodynamical approach to the nuclear dynamics in a 
superfluid fermionic environment 

2.1. Single nucleus in a superfluid medium 

Let us consider a spherical nucleus of a constant density pin and of 
radius R in the fluid at the density pout- The nuclear matter, both inside 
and outside a nucleus is assumed to be incompressible and irrotational, 
which implies that there exist a velocity potential <I>(r) = ^in{r) + $out(^) 
where $m(^ and $owt(r) stand for the velocity field inside and outside a 
nucleus, respectively. As a consequence of the continuity requirement, the 
function $ fulfills the Poisson equation: 

V'^^out = 0; for r > R, (2) 

V2$i„ = 0; for r <R. (3) 

The velocity field v{f) = Vin{r) + Voutir) is related to the velocity po- 
tential through the well known equations: 

Vin = V$in; for r> R, (4) 

Vout = "^^out; for r < R. (5) 
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The boundary conditions are determined by the geometry of the system. 
In the case of a spherical nucleus moving with the velocity u the following 
conditions have to be imposed on <I>: 



9^ 

— 5>i„ -n-u 
or 



^in\r=R — ^out\r=R, 

(6) 



r=R 

^out\r—KX) — 0, 



Pout ( '^^out - n-u 



r=R 



The vector n is the outward normal to the surface of a spherical nucleus 
and r is a radial coordinate (the center of a nucleus is placed at r = 0). 
The first equation stems from the requirement that the phase of the wave 
function of the superfluid system is continuous, the second one expresses a 
conservation of mass in radial flow, and finally the third one ensures that 
the correct asymptotic behavior for the velocity field is obtained. 

The solution of eq. (HJ with the boundary conditions © leads to the 
following expressions for the velocity potential: 

1-7 

^in(r) = -u-f (7) 

^ ^ 27 + 1 ^ ^ 

where 7 = Note that in the limit of 7 ^ 00, which correspond to the 

Pin 

case of an empty bubble moving through the fluid, one recovers the well 
known result PHIIHO], 

Thus $ describes a dipolar flow of the liquid around the impurity. The 
kinetic energy of the liquid associated with the flow reads: 



-m [ plV^Pd^r 
2 Jv ^ ' 



2 Jv 

1 /rfi ^^^"^Q 1 /rft 9^ out ^ a .Q^ 

= -"^Pm 9 S>i„— — - -mpout f 5>o„t— — di, (9) 

2 J s or 2 Js or 

where m is the nucleon mass and S is the surface of the nucleus. The 
last term appeared as a consequence of the Green's theorem. Applying the 
boundary conditions one obtains: 

T = ^m{pin - Pout) ^in\r=RU ■ ndS. (10) 

Hence the kinetic energy reads: 

2 „ d-^^a 

T = -TiR^mpin- —u^, (11) 

3 27 + 1 ' ^ ' 
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and it is clear that the effective mass of the bubble is equal to 

4 o (1-7)2 (1-7)2 

3 ^ 27 + 1 27 + 1 ' ^ ^ 

where M denotes the bare mass of the nucleus. 

Note that the effective mass depends on the nucleon density. It is in- 
structive to consider various limiting cases of the above formula, namely: 

j^ren ^ -mpoufT^R^ = -M for 7 ^ oo (bubble), (13) 
3 2 

j^ren _ g for 7 = 1 (uniform system), (14) 

j^ren _ —jyip^^-j^j^ = ]\/[ for 7 = (nuclcus in a vacuum). (15) 

3 

The above expressions show that the nuclear mass for a nucleus immersed 
in a superfluid neutron liquid is decreased as compared to the bare mass M 
(since pin > Pout in the case of a nucleus in a low density neutron fluid). 
Moreover, the renormalized mass depends on the density of the surrounding 
neutrons. This suggests then that dynamical properties of the Coulomb 
crystal will be strongly affected especially in the deeper layers of the inner 
crust. Indeed, the plasma frequency in the inner crust changes as compared 
to the value calculated for the case: pout = (ion in the vacuum). Namely, 

= \ , where uj'lf''^ corresponds to the plasma frequency for the 

system with renormalized ion masses |31| . 



2.2. Medium-induced interaction of two moving nuclei in a superfluid liquid 

Although the mass renormalization effect discussed in the previous sub- 
section seems to influence the dynamics of the ion crystal in the inner crust 
of neutron star, the picture is still not complete without taking into account 
the medium- induced interaction with other ions. Namely when the single 
nucleus is moving through the neutron liquid it produces a neutron flow, 
which affects the motion of neighbouring ions as well. One may expect that 
due to the fact that the potential flow generated by a moving nucleus be- 
haves like l/r^, it may induce a long range correlations in the system, and 
would supplement the standard Coulomb interaction between ions. 

Let us consider two nuclei moving with velocities Ui {i = 1,2). The 
nuclei have radii Ri and are placed at the distance r (see Fig. 1). The 
velocity potential can be expressed in the form: 

2 2 

$ = ^(^>i in + out) = ^^i in + ^out, (16) 
i=l 1=1 
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out 



i out J 



(17) 



i=l 



where $j in is the velocity potential inside the i—th nucleus and $j out is 
the outside potential depending only on the velocity Ui. The boundary 
conditions which have to be fulfilled in this case read: 



Pir. 



— $ 

dr ' 



in\r=Ri 

Hi ■ Ui 



r=Ri 
^out\i — >C» ~ 0) 



^out\r=Ri-i 

f 9 ^ 
Pout {-Q^^^out 



Hi ■ U, 



r=Ri 



(18) 



for i=l,2, 



where Ui is the outward normal to the surface of the i-th nucleus. We 
have also assumed here that the density inside each nucleus is the same and 
equal to pin- It is a justified assumption although it is easy to generalize 
the following expressions for the case of two diflFerent nuclear densities. 



Fig. 1. Two nuclei of radii and velocities and R2,U2, respectively 

The problem can be solved perturbatively, using the method of images, 
provided the distance between the nuclei is large. Namely, if r = \ fi — is 
large as compared to Ri and R2 then the natural small parameter appearing 

R R 

in the expansion of the velocity potential is — ^ — • Up to the second order 

in the expansion, the velocity potential outside nuclei can be expressed in 
the form: 

^ . fl • Ml 

^lOUt = Al g 
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+ 0[^A (19) 



r2 ■ U2 

^2 out = ^2 q — 



1/4 , X A \ f ^'^\ /^^1""2 3 (f- ■U2)(fl • r) 



+ -(^2 + ^2) — 



2 \ / \ ?"i f f 

+ (20) 

It is easy to check that, with the above choice of the outside velocity field, 
the following hmiting values of the potential and its radial derivatives are 
obtained: 

^ I . ^1 • ""1 . r-U2 

+ (^A2 + -5A2^ fRiY fRi-U2 3 {r.U2){Ri-f^' 



2 2 J \ r J \ Rf r2 Rf J 
+ o(^) (21) 



^ I . R2- U2 ^ . r- Ul 

'Pout \r2=R2 = ^2 ^3 + ^1^:3- 

+ (Iai + -6A,^ (R2\^ (R2-U1 3 (r-^xi)(i?2Tl' 



2 ' 2 ')\r J \ R^ r2 Rl j 
+ O(^) (22) 



-^22 3 {f ■ U2){Ri ■ r)^ 



M2 — 



+ O { ^) (23) 



9 , ^ . R2-U2 
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^ \ r ) \ R\ i?| y 

+ O f4V (24) 



where Ri = Ri — . Hence the inside solution can be written as: 



^lin = Ai—^ A2 



3 



Rf r 



2 2 J \ r J \ Rf r2 Rf 
+ O (1) , (25) 



^2m = A2—^+Ai 
^2 



+ 



3 
2 



•-2 

+ O ( — ) ; for n < i?i and rs < i?2- (26) 



The above form of the inside solutions ensure that the first pair of boundary 
conditions (|18() is automatically fulfilled. The other boundary conditions, 
originating from the conservation of the mass fiow, determine the unknown 
coefficients Ai,SAi, namely: 



,3 7-1 



Ai = Rf- , (27) 

'27 + 1' ^ ' 



,3 7-1 



SA, = -3i?f (^:^- (28) 

The total kinetic energy of the system consisting of two nuclei immersed in 
the superfiuid medium reads: 

T = \mprn( I |V«>i,„|Vri+ / |V«>2m|Vr2 

^ \Jri<Ri Jr2<R2 

+ Impout I |V$o«ipd^n, (29) 

2 JVout 

where Vout = V t^{R\ + ^2) volume ocuppied by the outside 

3 

neutrons. Analogously as in the case of a single nucleus one can transform 
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the above equation using the Green's theorem and the boundary conditions: 

T = \m{pin - Pout) i ^lin\r=Rx{ui ■ ni)dSl 



+ \^{Pin - Pout) / ^'2 
^ JS2 



in\r=R2 



{U2 ■ n2)dS2. 



(30) 



Hence finahy we get: 



1 



ren 



+ Anmpout 



1-7 
27 + 1 



^ ^ R1R2 



UI-U2 ^(ni • r)iu2 ■ r) 



where the renormahzed masses of nuclei have the form: 



Ml 



4 0(1 
3 ^ * 27 + 1 



Mi: 



27 + 1 ' 



(31) 



(32) 



for i = 1,2 and Mj denotes the nuclear bare mass of the i — th nucleus. 

As one can notice besides the usual kinetic energy of the two nuclei there 
appears also a crossterm, which in the case of two identical nuclei has the 
form: 



E.. 



int 




ui- U2 y{ui ■ r){u2 ■ r) 



UI-U2 7^{ui ■ r){u2 ■ r) 



(33) 



giving rise to the mutual interaction which is a manifestation of the Bernoulli 
law. The velocity dependent interaction between two nuclei has the same 
1/r'^ behavior as the (static) Fermionic Casimir interaction between two 
nuclei in a neutron Fermi sea, see Refs. |in[ IIH \T2 \ I13 | I14j. 

One may raise the following question: Does the medium-induced cou- 
pling between neighbouring nuclei H33|l represent a significant correction to 
the Coulomb interaction which stabilize the crystal? The problem is appar- 
ently innocuous since the interaction energy (|33|) depends on the velocity of 
nuclei, whereas the Coulomb energy is a position dependent function. Nev- 
ertheless, it is important to know whether the contributions coming from the 
hydrodynamic coupling represent only a small correction to the dynamical 
matrix, or may significantly modify the properties of the crystal. 

In the harmonic approximation one can estimate the ratio of the aver- 
aged hydrodynamic interaction energy and the Coulomb interaction energy. 
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associated with the thermal excitations 021- It turns out that the leading 

J R 

term is proportional to , where Rc is the radius of the Wigner- 

2j + lRc 

Seitz cell. Hence the hydrodynamical coupling turns out to be important 
in the bottom layers of the inner crust. 

3. Conclusion 

In the present paper we have formulated the hydrodynamic approach to 
the nuclear dynamics in the inner crust of neutron stars. We have derived 
the formula for the renormalized nuclear masses, assuming that the nuclei 
are immersed in an incompressible neutron superfluid. We argued that the 
renormalization of the mass affects the phonon spectrum of the crystal in 
a significant manner. Moreover, we have calculated the velocity dependent 
induced interaction between two moving nuclear impurities, which is a man- 
ifestation of the Bernoulli law. This interaction becomes important inside 
the inner crust. 

Summarizing, our investigations indicate that the dynamical properties 
of the nuclei in the neutron star inner crust are noticeably modified by 
the presence of superfluid neutron medium. One may expect that besides 
the translational degrees of freedom of the ions, also the nuclear shape 
vibrations can significantly influence the dynamic of the Coulomb crystal. 
Clearly, various nuclear collective degrees of freedom may turn out to be 
important for understanding of thermal and electric conductivities of the 
crust. 
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